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We analyze the ultra high resolution laser angle resolved photo-emission spectroscopy (ARPES) 
intensity from the slightly underdoped Bi2Sr2CaCu208+5 in the superconductive (SC) state. The 
momentum distribution curves (MDC) were fitted at each energy lj employing the SC Green's 
function along several cuts perpendicular to the Fermi surface with the tilt angle 9 with respect 
to the nodal cut. The clear observation of particle- hole mixing was utilized such that the complex 
self-energy as a function of to is directly obtained from the fitting. The obtained angle resolved 
self-energy is then used to deduce the Eliashberg function the diagonal channel by 

inverting the d-wave Eliashberg equation using the maximum entropy method. Besides a broad 
featureless spectrum up to the cutoff energy loc, the deduced a^F exhibits two peaks around 0.05 
eV and 0.015 eV. The former and the broad feature are already present in the normal state, while 
the latter emerges only below Tc. Both peaks become enhanced as T is lowered or the angle 9 moves 
away from the nodal direction. The implication of these findings are discussed. 

PACS numbers: 



I. INTRODUCTION 

The recent observation of the particle-hole mixing 
in the superconductive (SC) state of the cuprates 
by high resolution angle-resolved photo-emission spec- 
troscopy (ARPES) has opened up a new window to 
probe the fundamental physics of high temperature 
superconductivity.^'^ In particular, an analysis of the 
spectra in the SC state, using the Eliashberg formalism 
for d-wave superconductivity, provides the fluctuation 
spectrum responsible for pairing. This is an extension 
of the tunneling experiments and analysis with which it 
was definitively established that the pairing in metals like 
Pb is through exchange of phonons.^ It should be remem- 
bered that to get reliable information, it was necessary 
to have measurements of conductance at different tem- 
peratures and range of voltages of the order of the cut-ofF 
energy in the phonon spectrum to an accuracy of 0.2 %. 
The particle-hole mixing in cuprate superconductors was 
first observed some 15 years ago in ARPES<^ Those ex- 
periments had much worse momentum and energy reso- 
lutions. Since the cut-off is an order of magnitude higher 
for the cuprates than Pb and the angle-dependence of the 
spectra is crucial, the demands on the quality of the data 



are only being recently met through ultra-high resolution 
and stability of laser based ARPES. 

The ARPES provided an early evidence for the (1^2 _y2 
pairing state of the cuprates.— The measured leading edge 
shift of the energy distribution curve (EDC) of ARPES as 
a function of the tilt angle showed that the superconduct- 
ing gap is consistent with the d-wave pairing gap. It ush- 
ered in more debates and experiments which eventually 
led to the establishment of the d-wave pairing symme- 
try for the cuprate superconductors.^ The ARPES con- 
tains more information than the leading edge shift which 
may be utilized, for example, to extract the Eliashberg 
functions and track their evolution as the temperature 
is lowered below Tc. By properly extending the normal 
state analysis of extracting the self-energy, one should 
be able to deduce information about superconductivity 
of the cuprates such as the angle-resolved diagonal and 
off-diagonal self-energies, and the pertinent Eliashberg 
functions. This is precisely what we wish to present in 
this paper. 

For this, we fitted the ARPES momentum distribu- 
tion curves (MDC) at each energy uj employing the SC 
Green's function along several cuts perpendicular to the 
Fermi surface with the tilt angle 9 from the nodal cut 
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with respect to the (tt, tt) in the BriUouin zone. The 
clear observation of particle-hole mixing was utilized such 
that the complex self-energy as a function of lj is directly 
obtained from fitting the ARPES data. Thus obtained 
angle resolved diagonal self-energy Y,{9,uj) is then used 
to deduce the Eliashberg function a^F^^^^ (6*, w), i.e., the 
bosonic fluctuation spectrum multiplied by the coupling 
constant squared, in the diagonal channel by inverting 
the d-wave Eliashberg equation using the maximum en- 
tropy method (MEM). The diagonal self-energy evolves 
smoothly into the normal state self-energy as the temper- 
ature is raised above Tc- The evolution of the Eliashberg 
function as the temperature or tilt angle is varied will 
reveal a useful information about the nature of super- 
conductivity in the cuprates. 

On the other hand, the angle resolved off-diagonal self- 
energy (p{9,Lu), or, the density of states N{6,uj) given 
by Eq. ([6|) below, can be used for d-wave superconduc- 
tors to extract the Eliashberg function in the off-diagonal 
(i.e., pairing) channel^ as the ordinary tunneling conduc- 
tance was used by McMillan and Rowell to extract the 
spectrum of fluctuations for s-wave superconductors^ 
While the Eliashberg functions along the diagonal and 
off-diagonal channels are assumed to be the same for 
the s-wave pairing, they are in general different for d- 
wave superconductors. The current approach has the 
unique advantage in that it can disentangle the Eliash- 
berg functions in the diagonal and off-diagonal channels, 
a^F(+)(0,aj) and F^^'' {6 , uj) , respectively. 

In the following section II, we will present the formu- 
lation of the MDC analysis of the ARPES intensity in 
the superconducting state using the full momentum and 
energy dependence of SC Green's function. It is an ex- 
tension of the ARPES analysis in the normal state.— The 
results for the diagonal self-energy T,{9,uj) from slightly 
underdoped Bi2212 will be presented in section III at 
temperatures above and below Tc and along several cuts 
of the tilt angle 9 with respect to the (0, 0) — (tt, tt) nodal 
cut. As in the normal state, the extracted self-energy 
may be used as an input to deduce the Eliashberg func- 
tion a^F^+' (9, Lo). The obtained Eliashberg functions are 
presented in section IV. Recall that a^F(+^(6',w) along 
different cuts collapse onto a single curve with a peak 
near 0.05 eV below the angle dependent cutoff u!c{9) in 
the normal statCi^ In the SC state the peak around 0.05 
eV gets enhanced and additional peak emerges around 
0.015 eV below Tc. Both peaks become enhanced as T is 
lowered or the tilt angle is increased. We will conclude 
the paper by making some remarks and outlooks in the 
section V. 



where M(k, v) is the matrix element, v the energy of 
incident photon, /(w) the Fermi distribution function, 
A(k, oj) the quasiparticle (qp) spectral function, and 
B(k,uj) is the background from the scattering of the 
photo-electrons. We write the in-plane momentum k 
with the distance from the (tt, vr) point kj_ and the tilt 
angle measured from the nodal cut 9 as shown in Fig. 
[T] The self-energy has a much weaker dependence on k± 
than or a; as will be discussed below. Assuming this, 
the spectral function is written as 



A(k,uj) = ImG(k,w), 

TT 

W{9,u;)+Y{k,uj) 
W^{9,uj) -Y^{k,uj) - (j)^{9,uj)' 



(2) 



where (^(k, uj) is the retarded Green's function in the su- 
perconductive state. The following notations were used: 



W{9,u}) = u}Z{9,u}) = w - T.{9,Lo), 
Yik,Lo) = a^)+X{9,Lu), 
cl){9,io)^Z{9,oj)A{9,uj), 



(3) 



where Z{9,uj) is the renormalization function, X{9,ijj) 
the shift of the qp dispersion, and l](0,a;) and 4'{9,u)) 
represent the qp diagonal self-energy and the off-diagonal 
self-energy, respectively. The equation that connects 
S(0,w), X(0,w), and 0(6*, w) with the effective interac- 
tion in the charge and spin channels is the Eliashberg 
equation."^ It is presented in section IV below in connec- 
tion with extraction of the Eliashberg functions. 

It is informative to make the following decomposition 
of the SC Green's function: 



Y + W 



r2 _ 



l/2-|-A^/2 ^ 1/2- N 12 



Y 



Y + P 



(4) 



where 



P{9,uj) = ^/W'^{9,uj)-(t)^{9,u), 
W{9,Lo) 



N{9,uj) 



^W^{9,uj)-(t)^{9,uj)' 



(5) 
(6) 



II. FORMALISM 

The ARPES intensity, within the sudden approxima- 
tion, is given by 

I{k,uj) = \M{k,v)\'f{u:) [A{k,uo) + B{k,u:)] , (1) 



We note that the qp dispersion shift X{9,ijj) vanishes 
in the particle-hole symmetric band. Although the sym- 
metry does not hold for the realistic tight-binding disper- 
sion, it holds to a good degree over the small energy scale 
of SC and the renormalization is neglected in the present 
work. The ARPES intensity devided by the Fermi dis- 
tribution function is then given by 
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l + N{9,Lo) 



1 - N{0,LO) 



(7) 
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FIG. 1: The Fermi surface of Bi2212 in the Brillouin zone. 
The blue solid curves centered around the F point is the FS 
from Eq. (4HJ and the solid dots are the experimentally deter- 
mined FS at 6* = 0°, 5°, 10°, 15°, 20°, 25°. ifcx is the distance 
from the (vr, tt) point. The yellow thick curves along each cut 
indicate the actual momentum paths at w = of the experi- 
mentally measured ARPES MDC data. 



i{e,k±,uj) 

— = C{e,uj)lm 

.f w 



where C{d,uj) is the weight of the spectral function of 
the ARPES intensity. We then have the six parameter 
fit in SC state: C, B, the real and imaginary parts of P 
and as a function of binding energy w, while the nor- 
mal state fitting required four parameters. Note that the 
dependence on k± comes in through the bare dispersion 
^(k) only. It is therefore important to take appropriate 
dispersion. 

As in the normal state, we used the tight-binding (TB) 
dispersion and the linear dispersion (LD) for the MDC 
analysis. The TB dispersion ^(k) is given by 

^{kx, ky) — — 2i(cos kxa + cos kyu) + At' cos kxU cos kya 

— 2i"(cos 2kxa + cos 2kya) — /i,(8) 

where a = 3.82 A is the lattice constant and fi is the 
chemical potential. We took t = 0.395, t' = 0.084, 
t" = 0.042, and /i = -0.43 eV. The hnear dispersion 
was determined by linearization of the TB at FS of the 
six tilt angles 9. 

ao,k^)^VF{0)[k^~kF{0)], (9) 

where vp{6) and kp(9) are Fermi velocity and Fermi mo- 
mentum, respectively. The experimentally determined 
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FS in comparison with that from Eq. ([SJ is shown in Fig. 
[1] The six cuts with the tilt angles 6 with respect to the 
(tt, tt) are also shown with the solid lines. 

Note that the first and the second terms in Eq. ([7]) give 
the intensity due to the "particle" and the "hole" parts 
of the Bogoliubov particles, respectively. P{9,uj) and 
N{d,uj) are directly extracted from fitting the ARPES 
MDC data. Then, 

i:{0,ij) = uj-p{6,uj)N{e,u), (10) 

cj)\0,Lo) = P^{9,uj)[N\0,Lo)~l] (11) 

gives 1^(0, uj) and 4'{0,uj). Since the density of states 
N{0,uj) = 1 in the normal state, the off-diagonal self- 
energy (j){0, w) can only be extracted from the difference 
of the spectra between the normal and SC state. These 
differences are very small at energies above a few times 
Tc- So the requirements on the ARPES data to reli- 
ably extract 4> at higher energies are considerably more 
stringent than those to extract S. We defer this to future 
work and show here that considerable information on the 
fluctuation spectrum can be extracted from the diagonal 
self-energy alone. 
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FIG. 2: The representative MDC as a function of the mo- 
mentum along the tilt angle 6 — 20° . The dots are the exper- 
imental data and the solid red lines are the fitting. The first 
and second columns show the fittings in the normal state and 
in the SC state, respectively. The last column is the MDC 
ratios of SC to normal states. 



III. THE MDC ANALYSIS 

The ultra high resolution Laser ARPES data were 
collected from slightly underdoped Bi2212 of SC criti- 
cal temperature Tc — 89 K and pseudogap temperature 
T* « 160 K. The data were took along the cuts of the 
tilt angle 9 = (nodal cut), 5, 10, 15, 20, and 25 degrees 
with respect to the nodal direction and at temperatures 
T = 107, 97 above and 80, 70, and 16 K below T^. 
The photon energy of hiy = 6.994 eV was used in the 
laser ARPES. Refer to Ref.* for more detailed descrip- 
tion of the experimental setup and the technical details. 

Typical results of the ARPES analysis in SC state are 
given in Fig.[5]for the tilt angle 9 = 20°. The first column 
shows the results in the normal state at T = 97 K as 
a function of the magnitude of the inplane momentum 
k± from the (7r,7r) point at the energy w = —0.0005, — 
0.0205, and —0.1005 eV. The three energies represent the 
cases of w <C A, w sa A, and uj ^ A, where A is the gap 
amplitude at 61 = 20° and T = 16 K, 



A = 



Z(6l,A)' 



(12) 



The symbols are the data and the red solid lines are the 
fitting results. The agreements are almost perfect which 
justifies the neglect of the k± dependence of the self- 



energy. The second column is the corresponding results 
deep in the SC state at T = 16 K. The green solid (blue 
dashed) lines are the particle (hole) branch of the fitting, 
the first (second) term of Eq. ([7]). 

The important point is that the hole branch repre- 
sented by the blue curves exhibits a peak as a conse- 
quence of the particle-hole mixing of the pairing. This 
can be most spectacularly seen near w w A presented 
in the middle row. In addition to the main peak near 
k±a/TT « 0.91 from the original qp branch, there ex- 
ists the secondary peak at k±a/TT « 0.87. This is a 
direct observation of the particle-hole mixing deep in 
the SC state. The details of observations are presented 
separatelyji The particle-hole mixing was previously re- 
ported in the EDC by observing the bending-back of the 
spectral peaks."* Both branches of the Bogoliubov disper- 
sion due to the particle-hole mixing were also reported 
by the EDC in the intermediate temperature regime^ be- 
cause in the low temperature limit the Fermi function 
cuts the hole branch off and close to Tc the pairing fea- 
ture is very weak. The mixing is observed in the low 
temperature regime here and will be utilized to obtain in- 
formation about superconductivity in the cuprates. The 
last row shows the case of w ^ A. As the energy in- 
creases above A, the hole branch contribution vanishes 
as the bottom plots show. The last column is the ratios 
of the MDC at 16 K to 97 K which show the hole branch 
more clearly. 

We now show the real part of the extracted self-energy 
along the tilt angle 9 = and 9 = 20° in Fig. ^a) and 
(b), respectively, and that at T = 16 K in the plot (c). 
Fig. [3lja) and ^h) demonstrate that the feature around 
0.05 — 0.07 eV is already present in the normal state and is 
enhanced as the temperature is lowered, while the broad 
feature continues from the normal to SC state with no 
discernible change within the accuracy of the experiment. 
The — 0.02 eV feature emerges only along off-nodal cuts 
below Tc as can be seen from the plots (a) and (b), and 
its energy scale increases as the tilt angle is increased as 
can be seen from the plot (c). This is consistent with 
the d-wave pairing gap and implies that the — 0.02 eV 
feature is induced by superconductivity. 



IV. THE ELIASHBERG FUNCTION 

The extracted diagonal and off-diagonal self-energies 
may be used as experimental inputs to deduce the Eliash- 
berg functions F^'^^B , lo) and a^F^"' (6*, w) by invert- 
ing the Eliashberg equation. The d-wave Eliashberg 
equation may be written as 

oo />oo 

T.(9,uj)= I de de'S{uj,€,e')Ni{€)a'^F^+\9,e')(,13) 

OO J —OO 

OO pOO 

de / de'S{uj,e,e')Di{e)a^F^-\e')(,U) 



OO ^ — OO 



5(.,..') = M^,5) 

e + e' — u! — 10 
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focus on a^F(+'(0,a;) here. The real part of Eq. p3)) 
may be written as 




0.2 -0.1 

co(eV) 



FIG. 3: The real part of the extracted self-energy. The plots 
(a) and (b) are along the nodal cut and 6 — 20° , respectively, 
and the plot (c) is at temperature T = 16 K. 



where / and n represent the Fermi and Bose distribution 
functions, respectively. 
We took 

(l3{e,u) ^ (t){uj)sm{2e) (16) 
because the pairing is d-wave, and use the notations 

W{0',e) 



iVi(e) = Re 



1 



-Re 



(17) 
(18) 



^W^{e',e)-ct,^{e) sin2(20')/^, 
0(e) sin^ (2610 \ 
''^(^'^ ^Jw^^{e',e)~(j)■^{e) sin2(26")/^, 

a'F^^\e,e')^l'^^Fi+){e^e\e')) (,19) 
\ vf{9') / g, 



where vf{9') is the angle- dependent Fermi velocity and 
the bracket implies the angular average over 6'. 

As in the normal state, we invert the real part of the 
Eliashberg equation to deduce the Eliashberg function 
Q,2^(+) g'^ using the real part of the extracted diag- 
onal self-energy as an input. As mentioned before, the 
requirement on the data and numerical fitting are consid- 
erably stringent to reliably extract the off-diagonal self- 
energy than to extract the diagonal self-energy. We will 
defer the deduction of a^F^~\6,uj) to future work and 



Ei(0,cj) = 
K{uj,uj') 



dio'K{Lu,ijj')a^F^+\9,uj'), 
J{e)+n{-uj') 



de V- 



e + uj' — uj 



Ni{e) (20) 



where V represents the principal value, and the sub- 
scripts 1 and 2 refer to the real and imaginary parts. 
The inversion was performed using the maximum entropy 
methodj^ii^ Recall that in the normal state the Eliash- 
berg functions along different cuts of the tilt angle 9 all 
collapse onto a single curve which has a small peak at 
w 0.05 eV, flattens above 0.1 eV, and vanishes above the 
angle dependent cutoff uJc{0)- tOc{S) « 0.35 — 0.4 eV along 
the nodal direction and decreases as 6 increases 

Fig. 131 showing the deduced q;^F(+^ (6*, w), is the key 
results of the present paper. The deduced function is 
noisier at larger angles. The noise somewhat depends 
on the multiplier a of the maximum entropy methodJ^ 
We believe most of the oscillatory behavior seen are ar- 
tifacts of the MEM analysis and will focus only on the 
robust features whose variation is continuous as a func- 
tion of temperature and angles. The broad feature above 
about 0.07 eV does not change with angle or with tem- 
perature. It is a continuation of that required for the 
marginal Fermi liquid properties which were derived re- 
cently to arise from quantum criticality.— It carries about 
3/4 of the total spectral weight. 

The plot Ufa) is slightly below at T = 80 K. 
The Eliashberg function does not change much from 
the normal state shape except that the peak value at 
Lj « 0.05 eV increases to approximately 0.4 for large tilt 
angles from 0.3 of the normal state value. In Fig. Hljb), 
a'^ F^-^\9 , Lo) along the nodal cut is shown as the tem- 
perature is varied. As might be expected from above 
behavior, there is little change along the nodal cut, al- 
though there is a sign of the lower energy peak at T = 16 
K. It seems that the change in S(0 = 0,w) as T is varied 
as shown in Fig. [3l^a) is predominantly from the change 
in the density of states that enters Eq. (l20l) . 



Q,2^(+) (^9,uj) deep in the SC state is shown in Fig. ■ 
Like at T = 80 K, the peak at a; = 0.05 has its normal 
state value for small tilt angle and increases as the angle 
increases. Also a second peak at w 0.015 eV emerges 
which, like the w w 0.05 eV peak, increases from the 
normal state value as the angle is increased. In Fig.|4l^d), 
the Eliashberg function along the cut 6 = 20° is shown. 
As the temperature is lowered below Tc, the 0.05 eV peak 
is enhanced and the 0.015 eV peak newly develops, as one 
can anticipate from preceding discussion. Both peaks are 
enhanced as T is lowered or the tilt angle is increased. 

There have been many investigations of the fluctuation 
spectrum of the cuprate superconductors such as the in- 
frared conductivity, inelastic neutron scattering, Raman 
scattering, scanning tunneling spectroscopy, and so on, 
which are less direct in the information they provide for 
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FIG. 4: The deduced Eliashberg function a'^ F^+\e , uj) . Plot 
(a) is at T = 80 K slightly below Tc, plot (b) along the nodal 
cut, plot (c) at T = 16 K deep in the SC state, and (d) is 
along the off-nodal cut of 6 = 20°. 



the source of superconductivity than ARPES. Analysis 
of the frequency dependent conductivity by the McMas- 
ter group reported that a single peak shows up in the 
Eliashberg function below 0.1 eV for Bi2Sr2CaCu208+5 
and other cuprate compounds, but a double peak feature 
for Lai.83Sro.i7Cu04 at low temperatures. It exhibits a 
peak at a; « 0.05 eV at high temperature T — 250 K, 
but as T is lowered to 30 K it showed two peaks at 0.015 
and 0.044 eVJ^ The inelastic neutron scattering (INS) 
experiment on Lai.84Sro.i6Cu04 also reported the two 
peak structure around 0.018 and 0.04 — 0.07 eV at the 
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FIG. 5: The 9' averaged susceptibility Xii^y^) calculated 
from VignoUe to compare with a^F'~'^\6,uj). The plot (a) 
is the Xii^y^) from VignoUe and (b) is with the correlation 
length reduced to 0.1 of (a). 



antiferromagnetic wave vector. ^■^ The Eliashberg analysis 
of the break junction SIS conductance and scanning tun- 
neling spectra on overdoped Bi2Sr2CaCu208+(5 reported 
a single peak near « 0.02 eVM' 

No inelastic neutron scattering results are available 
for Bi2Sr2CaCu208+5, so that we can only compare 
our results with the detailed extraction of the spec- 
tral function of the magnetic fluctuations i^) for 
Lai.84Sro.i6Cu04. The positions of the two peaks around 
0.018 and 0.04 — 0.07 eV are consistent with the peak po- 
sitions of the deduced {9 , uj) . To make comparison 
of the momentum dependence between our a^F^^^OjUi) 
and INS results of X2(q, we compute X2{d7 ^) by tak- 
ing the integral over 6' with both k and k' on the Fermi 
surface having the tilt angles 6 and 9', respectively. As 
with Eq. (ITOl) . we take 



X2(^,c^) = (X2(k-k',w)), 



(21) 



The imaginary part of the spin susceptibility X2(q, 
was taken from VignoUe et al. INS results»i2, Fig. [Sfa) is 
the computed X2(^, ^) from VignoUe and (b) is that with 
the correlation length reduced to 1/10 of (a). Comparing 
with our results in Fig. 21 the variation with angle com- 
pares better with the reduced correlation length. The 
angle independence of the Uat part of F''~^\9 , lo) is 
consistent only with the fluctuation spectrum which has 
the short correlation length on the scale of the lattice 
constant. 

The origin of the ~ 50 meV feature is often taken to 
be a phononi^ii^ but it may be the strong dispersionless 
magnetic feature recently observed^^ii^ only in the pseu- 
dogap region in the cuprates. One may be able to decide 
between the two through ARPES analysis similar to that 
done here of overdoped samples without the pseudogap. 
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V. CONCLUDING REMARKS 

We have presented the analysis of the ARPES inten- 
sity in the superconductive state for the Eliashbcrg func- 
tion along the diagonal channel, F^~^\0 , u) . Beside 
the broad featureless spectrum, a peak at 0.05 eV present 
in the normal state is enhanced as T is lowered, and a 
second peak emerges around 0.015 eV in the SC state. 
The 0.015 eV peak is an interesting new feature which 
needs further exploration. Since it appears only below 
Tc, it can not be responsible for Tc. We do not have an 
answer at present to the question as to what kind of col- 
lective modes of the superconducting state it represents. 



The origin of the 50 mcV feature generally observed is 
at the heart of current debate. One may be able to an- 
swer through ARPES analysis similar to that done here 
of overdoped samples without the pseudogap. 
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